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INTERPOLATION  AND  SMOOTHING 
OF  EXPERIMENTAL  DATA 
WITH  SLIDING  POLYNOMIALS 

By  W.  M.  Snyder' 

ABSTRACT 

This  paper  presents  a  comprehensive  overview  of  the  analytical  structure 
of  sliding  polynomials.  Emphasis  is  placed  on  explanation  so  that  researchers 
with  limited  background  in  numerical  analysis  can  apply  the  methods. 
Presentation  is  generalized  to  assist  readers  from  various  disciplines  in 
preparing  data  for  interpretation.  Polynomials  of  increasing  complexity  were 
previously  presented  in  separate  publications.  The  intent  here  is  to  develop 
an  appreciation  of  the  common  elements  underlying  three  levels  of  complexi- 
ty, starting  with  simple  linear  forms.  Continuity  of  the  first  derivative  of 
interpolation  arcs  is  added  first,  followed  by  extension  to  continuity  of  the 
second  derivative.  The  three  forms  are  then  transformed  to  smoothing 
operators  using  least  squares. 


INTRODUCTION 

Interpolation  and  smoothing  are  routinely  re- 
quired to  reduce  and  process  experimental  data. 
Choice  of  a  particular  method  depends  on  the 
amount  of  data  and  on  awareness  of  available 
methods.  DuChateau  et  al.  (1972)  and  Kimball 
(1974)  dealt  with  interpolation  and  smoothing 
techniques  that  require  substantial  amounts  of 
experimental  data.  Such  data  may  be  generated 
either  in  an  extensive  cyclic  time  domain  or  by 
repetition  in  a  limited  time  or  space  domain. 
Savitzky  and  Golay  (1964)  discussed  smoothing 
and  differentiation  of  data  using  convoluting 
functions  in  weighted  moving  average  tech- 
niques. Again,  availability  of  substantial  a- 
mounts  of  data  is  implied,  generated  by  basically 
continuous  physical  experimentation.  The  author 
(1961,  1962,  1967a)  showed  how  sliding  polyno- 
mials could  be  used  for  interpolation  and  smooth- 
ing. This  technique  can  be  used  with  more 
limited  amounts  of  experimental  data. 

1  Hydraulic  engineer.  Southeast  Watershed  Research 
Laboratory,  Agricultural  Research  Service,  U.S.  Depart- 
ment of  Agriculture,  Athens,  Ga.  30601. 


This  publication  draws  together  and  expands 
the  author's  previous  work  on  sliding  polyno- 
mials. The  construction  of  sliding  polynomials  is 
reviewed  and  analyzed.  Common  elements  and 
concepts  can  be  found  in  a  hierarchy  of  polyno- 
mials. Tracing  these  common  elements  through 
three  levels  of  complexity  will  be  useful  to 
researchers  in  many  disciplines  who  must  pro- 
cess observational  or  experimental  data.  Realiza- 
tion of  the  full  potential  of  the  sliding  polynomials 
will  aid  in  choosing  interpolation  and  smoothing 
techniques  best  suited  to  specific  data  processing 
requirements. 

Understanding  the  hierarchial  structure  of  the 
polynomials  can  begin  with  review  of  simple 
linear  interpolation.  Extension  of  this  simple 
form  to  piecewise  linear  functions  establishes  a 
basis  for  understanding  more  complex  forms. 
Addition  of  a  requirement  for  mathematical 
continuity  of  a  first  derivative  is  a  logical  next 
step.  Further  extension  follows  with  require- 
ment for  mathematical  continuity  of  a  second 
derivative.  All  of  these  three  piecewise  forms  are 
utilitarian.  While  a  basic  knowledge  of  mathe- 
matical and  numerical  analysis  is  assumed,  em- 
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phasis  is  placed  on  explanation  to  aid  the  reader 
in  selection  of  the  form  suited  to  his  purposes. 

INTERPOLATION 
WITH  POLYNOMIALS 

Interpolation  is  the  estimation  of  values  be- 
tween given  discrete  values  of  a  property.  Some 
function  must  be  assumed.  This  function  must 
give  values  of  the  property  dependent  upon 
values  and  locations  of  the  given  discrete  points. 
A  one-to-one  correspondence  exists  between 
mathematical  degrees  of  freedom  in  the  function 
and  the  number  of  given  discrete  points.  Data 
points  through  which  the  function  is  required  to 
pass  will  be  called  base  points. 

In  the  following  discussion,  basic  polynomial 
properties  and  the  functional  forms  of  sliding 
polynomials  are  explained.  Then,  the  use  of  these 
concepts  in  interpolation  is  described. 

Basic  Polynomial  Properties 

Polynomials  of  any  degree  could  be  assumed 
for  interpolation.  Beyond  the  fifth  power,  how- 
ever, the  algebraic  manipulations  become  tedious 
and  error  prone.  A  fifth-power  polynomial, 
evaluated  from  six  or  more  coordinate  points,  is 
sufficient  to  produce  smoothly  continuous  rates 
of  change  through  experimental  data.  Therefore, 
this  presentation  does  not  go  beyond  fifth -power 
polynomials. 

The  fifth  -power  polynomial  has  the  conven- 
tional form 

y=a+bx+cx2+dx3+ex4+fx5.  (1) 

Coefficients  a  through  /  represent  six  mathe- 
matical degrees  of  freedom  which  can  be  satisfied 
by  reference  to  six  given  discrete  points.  These 
concepts  are  developed  under  "Six-Point  Inter- 
polation." 

The  first  derivative  of  equation  1,  the  rate  of 
change  of  y,  is 


dy/dx=b+2cx+3dx2+'lex:i+5fxi  (2) 

Higher  order  derivatives  are  given  in  equa- 
tions 3  through  6. 

d2y/dx2=2c+6dx+\2ex  +20fx3.  (3) 

d?y/dx3=6d+24ex+60fx2.  (4) 

d4y/dx4=2ie  +  120fx.  (5) 

d5y/dx$=l20f.  (6) 


Note  the  fundamental  continuity  of  form  in 
equations  1  through  5.  Any  one  of  these  can  be 
used  as  an  interpolation  function.  Properties  oi 
the  function  for  estimating  rates  of  change  are 
given  by  derivatives  of  order  higher  than  tht 
chosen  functional  form.  Properties  of  the  func 
tion  for  integration  are  given  by  derivatives  of 
order  lower  than  the  chosen  functional  form.  For 
example,  equation  5,  the  straight  line,  has  a 
constant  rate  of  change  given  by  equation  6  and  a; 
parabolic  rate  of  accumulation  given  by  equation 
4,  In  the  discussions  of  two-point,  four-point,  and 
six-point  interpolation,  equations  5,  3,  and  1  are! 
used  to  develop  special  polynomials.  They  act  as  j 
flexible  splines  that  pass  through  extended  num- 
bers of  coordinate  points.  Since  symmetry  is 
required,  only  polynomials  with  even  degrees  of 
freedom  can  be  used  in  this  method. 


Sliding  Polynomials 

Different  terminology  may  be  used  to  describe 
the  mathematical  techniques  of  passing  a  contin-  j 
uously  transhaping  function  from  left  to  right 
through  some  set  of  discrete  values,  or  base 
points.  As  the  function  moves  to  the  right,  base 
points  must  drop  out  on  the  left  side  of  the 
function  span,  and  new  points  must  be  picked  up 
on  the  right  side  of  the  function  span.  This 
concept  can  be  called  continuous  transforming, 
flexible  spline,  or  sliding  polynomial.  The  term 
"sliding  polynomial"  is  considered  most  suitable 
in  the  development  that  follows. 

First,  consider  the  schematics  in  figure  1.  They 
are  sufficient  for  a  simple  situation  where  inter- 
polation is  always  within  one  interval  between 
base  points. 

Note  the  set  of  coordinate  points  YBO  through 
YB6  in  figure  la.  These  coordinate  points,  spaced 
at  uniform  intervals,  will  be  used  as  base  points 
through  which  interpolation  polynomials  are  re- 
quired to  pass.  The  designation  YBO,  for  exam- 
ple, means  a  base  point  in  y  at  a  value  of  x  equal 
to  zero.  Uniform  spacing  of  the  base  points  is  not 
absolutely  necessary,  but  it  does  greatly  simplify 
the  algebraic  forms. 

Now  observe  in  detail  the  interval  between 
YB2  and  YBS.  Figure  lb  is  the  schematic 
representation  for  passing  a  straight  line  be- 
tween the  two  points.  Interpolation  can  be  made 
on  this  straight  line.  Figure  lc  shows  a  cubic 
equation  passing  through  YB1,  YB2,  YBS,  and 


2 


YBO 
O 


c. 


YBI 
O 


YB2 
O 


YB3 

O 


YB5 

O 


YB6 

O 


YB4 

O 


Scale  of  X 


0  12  3  4  5 

a.    DEFINITION  OF  BASE  POINTS 


YB2 
O— 


YB3 
— O 


y  =  a  +  bx 
TWO- POINT  INTERPOLATION 


YB2_ 


YB3 


YBI  - 
O- 


^YB4 

3 


y  =  a  +  bx  +  cx2  +  dx3 
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y  =  a  +  bx  +  cx2  +  dx3  +  ex4+fx5 
d.    SIX-POINT  INTERPOLATION 

FIGURE  1 . — Definition  of  interpolation  systems. 


on  the  left  and  picking  up  YB5  on  the  right.  The 
phasing  out  and  phasing  in  of  the  base  points  is 
gradual,  instead  of  abrupt,  as  in  the  linear 
segments. 

Figure  2c  is  an  extension  of  concepts  in  figure 
2b.  The  higher  order  polynomial  phases  out  YBO 
and  phases  in  YB6  as  the  ar-point  leaves  the 
YB2-YB3  arc  and  passes  through  x=3  into  the 
YB3-YBA  arc.  Special  forms  for  this  method  are 
shown  under  "Six-Point  Interpolation." 

Two- Point  Interpolation 

Two-point  interpolation  is  the  familiar  and 
conventional  linear  interpolation.  The  method  is 
intuitively  understood  and  used  by  many  per- 
sons. This  intuitive  understanding  is  used  here  as 
the  basis  for  presentation  of  more  complex 
polynominal  methods.  To  prepare  for  this  pres- 
entation, the  familiar  linear  interpolation  is  cast 
into  algebraic  forms  which  will  be  needed  later. 
Only  two  base  points  at  a  time  are  used  in  linear 
interpolation.  In  figure  la,  the  two  points  are 
YB2  and  YB3.  The  two  parameters  in  the 
associated  equation  can  be  evaluated  in  terms  of 
the  base  points.  The  matching  of  number  of 
parameters  and  number  of  base  points  should  be 
noted. 


VBA.  Interpolation  is  to  be  limited  to  the  interval 
)etween  YB2,  and  YB3,  the  same  interval  cov- 
ered by  the  straight-line  example.  Figure  Id 
shows  a  fifth -power  equation  passing  through  six 
joints.  Interpolation  again  is  to  be  limited  to  the 
nterval  between  YB2  and  YB3. 

A  more  realistic,  and  more  complicated,  prob- 
em  is  interpolating  in  adjacent  intervals  (fig.  2). 
figure  2a  shows  interpolation  on  two  adjacent 
inear  segments.  Visualize  a  point  of  interpola- 
tion on  the  x-scale  moving  from  YB2  to  YB3. 
Values  of  y  can  be  read  from  the  straight  line 
;onnecting  YB2  and  YB3.  When  the  point  passes 
through  YB3,  a  transformation  is  necessary. 
Point  YB2  is  dropped  on  the  left,  and  point  YB4  is 
idded  on  the  right.  As  the  x-point  moves  through 
c=3,  interpolation  changes  abruptly  from  the 
ine  through  YB2  and  YB3  to  the  line  through 
YB3  and  FB4. 

A  smoother  transition  through  YB3,  utilizing 
in  arc  through  YB2  and  YB3  and  one  through 
YB3  and  YB4,  is  shown  in  figure  2b.  The 
:echnique  to  suppress  the  angle  between  linear 
segments  at  YB3  is  developed  under  "Four-Point 
nterpolation."  It  is  based  on  dropping  out  YBI 
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FIGURE  2. — Definitions  for  connected  arcs. 
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Define  a  new  scale  variable  z,  with  2=0  at 
x=2,  and  2=1  at  x=3.  An  equation  to  define  the 
assumed  linear  interpolation  function  between 
YB2  and  YBS  can  now  be  written: 

y=a+bz.  (7) 

At  2=0,  y=YB2,  and  2=1,  y=YB3.  Substi- 
tuting these  coordinate  values  of  the  two  base 
points  in  equation  7  yields  the  simultaneous 
equations 

YB2=a  (8) 
YB3=a+b.  (9) 

From  equations  8  and  9,  a=YB2,  and  b  =  YBS 
-YB2.  Equation  7  may  now  be  written  in  terms  of 
the  base  points  as 

y=YB2+(YB3-YB2)z=(l-z)YB2+zYB3.  (10) 

The  solution  for  equations  8  and  9  may  also  be 
written  in  the  form  of  the  matrix-vector  equation 
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So-called  interpolation  coefficients  can  be  de- 
veloped from  special  values  of  YB2  and  YBS. 
First,  set  YB2  equal  to  1  and  YBS  equal  to  zero. 
Expansion  of  equation  11  gives  a=l  and  6  =  -l. 
The  special  form  of  equation  10  results  in 

ya-\-z.  (12) 

Since  YB2  was  set  to  unity,  ya  can  be  called  the 
unit  contribution  of  a  base  point  at  2=0,  the  left 
one  of  two  base  points.  Next,  set  YB2  equal  to 
zero  and  YBS  equal  to  1.  Expansion  of  equation 
11  gives  a=0  and  6  =  1.  The  special  form  of 
equation  10  now  is 

yh=z.  (13) 

Now,  yb  can  be  called  the  unit  contribution  of  a 
base  point  at  2=1,  the  right  one  of  the  two  base 
points. 

Equations  12  and  13  can  be  evaluated  for 
values  of  z  as  in  table  1.  The  values  of  ya  and  yb 
are  normally  designated  as  CO  and  CI,  meaning 
the  interpolating  coefficient  to  be  multiplied  by 
the  base  value  at  2=0  and  by  the  base  value  at 
2=1.  The  same  meaning  is  expressed  on  the 
right  in  equation  10.  For  example,  two-point 
interpolation  at  a  spot  one-fourth  of  the  x-dis- 
tance  from  YB2  to  YBS  is  written  directly,  with 
coefficients  from  table  1,  as 

2/=(0.75)FB2+(0.25)yfi3.  (14) 

Equation  14,  of  course,  is  the  familiar  and 


conventional  form  for  linear  interpolation. 

If  the  two  adjacent  linear  segments  in  figure 
are  considered,  equations  can  be  written  direc 
for  interpolating  y  x  on  the  left  segment  and  y2 
the  right  segment: 


n 


and 


y  =YB2+(YB3-YB2)z, 
y  =  YB3+(YB4-YB3)z. 
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With  2=1,  equation  15  gives  ?/,=FZ?3. 
with  2=0,  equation  16  gives  y2  =  YBS.  T 
simply  confirms  that  both  segments  pass  throuki 
YBS,  and  the  interpolation  is  continuous  in  t» 
sense  that  there  is  no  gap  at  the  base  point  Ylh 
However,  dy  J dz=YBS -YB2  and  dy2/dz=Y% 
-YBS  show  that  the  slopes  of  the  segments  a»|) 
different,  constant,  and  not  continuous. 


TABLE  1  —  Two-point  interpolation  coefficient  \ 


z 

Va 

Vb 

0.10 

0.90 

0.10 

.25 

.75 

.25 

.50 

.50 

.50 

.75 

.25 

.75 

.90 

.10 

.90 

Four-Point  Interpolation 

In  geometrical  concept,  two-point  interpol 
tion  can  be  visualized  as  a  point  being  linear1 
displaced  in  t/-scale  as  it  moves  along  the  2-sca 
from  one  base  point  to  the  next.  This  tw 
directional  displacement  generates  a  locus 
interpolated  values  which  lie  along  a  straigi 
line.  Note  particularly  the  concept  of  a  point 
continuous  displacement  generating  a  straigh! 
line  interpolation  locus,  and  note  that  the  line 
one  mathematical  order  up  from  the  point. 

In  four-point  interpolation,  the  geometric 
concept  is  a  function  that  will  be  linearly  tra; 
shaped  in  ?/-scale  as  a  point  of  interpolate 
moves  along  the  2-scale.  This  function  is  up  tv 
orders  from  the  point,  and  is,  therefore, 
parabola.  The  locus  of  interpolated  values  is  oi 
order  higher  than  the  function  being  transhap( 
and  is  a  cubic  polynomial. 

The  geometrical  concept  on  which  four-poii 
interpolation  is  based  is  further  developed  ; 
follows:  Visualize  a  parabola  that  passes  throug 
base  points  YBl,  YB2,  and  YBS  in  figure  1.  Nex 
visualize  the  parabola  through  YB2,  YBS,  ar 
YB4.  As  the  point  of  interpolation  moves  alor 
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llhe  x-scale  from  YB2  to  YB3,  the  interpolated 
italue  is  to  be  read  from  a  flexible  parabola  that 
Kranshapes  linearly  from  the  left  parabola  to  the 

ight  parabola. 
The  analytical  concept  of  four-point  interpola- 
tion requires  specification  of  four  mathematical 
:  onstraints  on  which  solution  for  the  four  coeffi- 
cients of  the  cubic  polynomial,  equation  3,  are  to 
he  based.  The  mathematical  specifications  are 
ijased  on  commonality.  The  cubic  equation  is 
required  to  be  common  to  the  parabolas  at  the 
imds  of  the  interpolation  arc.  This  means  simply 
lihat  both  the  parabolas  and  the  cubic  pass 
jhrough  YB2  and  YB3  (fig.  lc).  In  addition,  the 

parabolas  and  the  cubic  are  required  to  have  the 
fame  slope,  or  the  same  first  derivative,  at  YB2 

,nd  YBS. 

I  Next,  specify  a  scale  variate  2,  which  takes 
falues  0,  1,  2,  and  3  from  left  to  right  in  any 
pur-point  set  of  base  points.  Define  the  left 
arabola  by 

j  y=u  +  VZ+WZ2  (17) 

|ad  d1y/dz=v+2wz.  (18) 

\  Solve  for  coefficients  u,  v,  and  w  by  passing 
irough  the  three  points  2=0,  1y=YBl;  2=1, 
y=YB2;  and  z=2,  Xy—YB3,  giving  the  values 

u=YBl,  (19) 
l  v  =-(3/2)YBl+2-YB2-<l/2)FB3,  (20) 

id  u,=(l/2)YBl-YB2+(l/2)FB3.  (21) 

1  Next  define  the  right  parabola  by 

2y=u+vz+wz  (22) 
,nd  d2y/dz=v+2wz.  (23) 

|  Solve  for  new  coefficients  u,  v,  and  w  by 
(assing  the  parabola  through  points  2  =  1, 
y=YB2;  z=2,  2y=YB3;  and  2=3,  2y=YB4. 

'he  solution  is 


u=3-YB2-3  YB3+YB4,  (24) 

u=-(5/2)VB2+4VB3-(3/2)YB4,  (25) 

id             u;=(l/2)YB2-YB3+(l/2)YB4.  (26) 

Define  the  cubic  equation  by 

y=a+bz+cz  +dz3  (27) 

id                   dy/dz=b+2cz+Mz2.  (28) 


The  four  mathematical  constraints  expressing 
lommonality  between  parabolas  and  cubic  can 
low  be  stated  as  four  simultaneous  equations 
29-32).  Equations  17  and  18,  with  coefficients  in 
equations  19,  20,  and  21,  are  set  equal  to 


equations  27  and  28  for  2=1.  Equations  22  and 
23,  with  coefficients  in  equations  24,  25,  and  26, 
are  set  equal  to  equations  27  and  28  for  2=2.  The 
resulting  four  simultaneous  equations  are 


YB2=a+b+c+d,  (29) 

(l/2)yB3-(l/2)YBl=&+2c+3eZ,  (30) 

YB3=a+26+4c+8d,  (31) 

and  (l/2)YB4-(l/2)YB2=c-+4c+12d.  (32) 

The  solution  to  these  equations  is 

a=2-YBl-3-YB2+3-YB3-YB4,  (33) 

&=_4.yBl+(19/2)yB2-8-YB3+(5/2)yB4,  (34) 

c=(5/2)FSl-7-yB2+(13/2)yB3-2-VB4,  (35) 


and    d=-<l/2)yBl+(3/2)YB2-(3/2)VB3+(l/2)yB4.  (36) 

This  solution  can  be  expressed  concisely  in  the 
matrix-vector  equation 
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Next,  consider  the  property  of  the  locus  of 
interpolated  values  at  the  common  point  of  two 
arcs.  Specifically,  in  figure  2b,  arcs  YB2-YB3  and 
YB3-YBA  obviously  join  at  YB3.  Now  at  YBS, 
interpolation  arc  YB2-YB3  ends  on  the  parabola 
through  YB2,  YB3,  and  FS4.  However,  interpo- 
lation arc  YB3-YB4  takes  off  from  this  same 
parabola.  The  four  mathematical  specifications 
require  the  interpolation  cubics  to  be  common 
and  tangent  to  the  parabolas  at  the  ends  of  the 
arcs.  The  cubics  of  the  YB1-YB4  set  and  the 
YB2-YB5  set  are,  therefore,  both  common  and 
tangent  to  the  same  parabola  at  YB3.  The  cubics 
of  adjacent  interpolation  arcs  are  thus  common 
and  tangent  to  each  other. 

Thus,  it  has  been  demonstrated  that  four-point 
interpolation  based  on  a  sliding,  transhaping 
parabola  generates  values  on  a  cubic  polynomial. 
Adjacent  cubic  arcs  are  continuous  and  have 
continuous  first  derivatives  as  the  interpolation 
locus  passes  through  the  base  points. 

Figure  3  is  an  example  of  four-point  interpola- 
tion. (Computer  program  A-l,  appendix  A,  gen- 
erated the  interpolated  values.)  Remember, 
when  looking  at  figure  3,  that  each  interpolation 
arc  is  a  different  cubic  polynomial,  with  the 
ensemble  of  cubics  having  the  unique  property 
that  they  are  common  and  tangent  at  the  base 
points. 
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FIGURE  3. — Four-point  interpolation  and  integration. 

In  summary,  four-point  interpolation  estab- 
lishes a  form-free  ensemble  with  two  orders  of 
continuity.  Generally,  the  interpolation  locus  will 
pass  through  the  base  points  without  undesirable 
oscillation.  In  some  situations,  with  an  abrupt 
change  in  slope  of  the  base  point  data,  a  tendency 
to  oscillate  may  develop. 

Four-point  interpolation  can  be  reduced  to 
coefficients  that  operate  directly  on  the  base 
points,  just  as  the  two-point  system  was  re- 
duced. Consider  the  special  case  of  YB\  —  \  and 
YB2=YB3=YB4=0.  Then,  from  equations  33-36 
or  equation  37,  a=2,  b=~4,  c=5/2,  and  d— 
_l/2,  and  the  interpolation  cubic  is 

j/0=2-4r+(5/2)22-(l/2)23.  (38) 

Equation  38  can  be  evaluated  for  selected 
values  of  z,  producing  a  coefficient  that  gives  the 
unit  contribution  of  the  base  point  YB1  to  the 
interpolated  value.  Since  YB1  is  at  a  z-value  of 
zero,  call  this  coefficient  CO. 

Now  use  the  special  case  YB2=\  and  YB1 
=  YB3=YB4=0.  The  interpolation  cubic  from 
equations  33-36  is 

yi=3+(19/2)z-7z^+(3/2)z:i.  (39) 

Evaluate  this  for  the  same  selected  z-values 
used  with  equation  38.  Call  this  set  of  coefficients 


CI,  since  YB2  is  at  a  z-value  of  1. 

Similarly,  evaluate  C2  and  C3.  An  interpolate^ 
value  is  the  summation  of  the  products  of  the  uni  1 
operators  with  any  real  set  of  YBl  through  YB4{ 
This  product  summation  is 

y=yo  +  yi+y2+y.3-C0-YBl  +  ClYB2  +  C2-YB3 

+C3-YB4.  (40) 

Table  A-l  (appendix  A)  is  a  listing  of  the  foui 
interpolating  coefficients  for  each  one-hundredtl* 
of  the  interpolating  arc  from  2=1  to  z=2.  These! 
coefficients  were  generated  by  program  A-2. 

The  area  under  any  function  is  given  by  the 
integral  of  the  function  across  a  range  of  x  (or  z). 
The  integral  of  the  interpolation  cubic  of  the 
four-point  method  is  given  in  equation  41.  Inte- 
gration, just  as  interpolation,  must  be  limited  tc 
the  central  interval  of  the  four-point  set: 

V  ydz=az  +  (b/2)z2  +  {c/3)z3  +  (d/4)zA  +  C=A.  (41)! 

The  constant  of  integration,  C,  can  be  evalu- 
ated by  noting  that  the  area  under  the  curve,  M 
is  zero  when  z=l;  that  is,  at  the  left  end  of  the 
interpolation  arc.  For  this  value  of  z,  and  with 
equations  33-36,  C  is  given  by 

C=(l/24)(-17-YBl+5-YB2-19-YB3+7-VB4).  (42; 

The  area  under  any  segment  of  the  interpola- 
tion arc  can  now  be  calculated  with  equations 
33-36,  41,  and  42.  It  is  also  possible  to  reduce: 
four-point  integration  to  a  coefficient-multiplier 
scheme  analagous  to  equation  40  for  interpola- 
tions. Table  A-2  is  a  listing  of  such  coefficients, 
also  generated  by  program  A-2  (appendix  A).jj 
Now,  note  that  four-point  integration  is  also 
included  in  program  A-l.  The  area  under  the  arcs 
of  the  interpolation  cubic  are  plotted  in  figure  3. 
Remember  that  each  arc  starts  with  zero  area. 
The  total  area  to  the  left  of  any  point  is  thus  the 
sum  of  the  partial  area  in  the  arc  plus  all  full  arcs 
to  the  left. 

The  coefficients  in  table  A-2  are  based  on  a 
scale  variate  z  with  the  base  points  at  unit 
intervals  of  z.  If  the  base  points  are  at  a  scale  s  in 
the  original  data,  given  by  s=kz+z0,  then 
ds=kdz.  An  integral  ydz  can  be  changed  to  an 
integral  yds=kydz  by  simple  change  of  variable. 
This  change  of  variable  is  included  in  program 
A-l,  but  must  be  performed  as  a  separate 
operation  following  use  of  table  A-2. 
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Six-Point  Interpolation 

The  development  of  the  method  of  six-point 
interpolation  has  been  given  in  detail  previously 
[(Snyder  1967a).  The  geometrical  and  analytical 
[concepts  are  analagous  to  those  used  in  develop- 
ment of  the  four-point  method  in  the  previous 
section. 

[  Note  briefly  that  six-point  interpolation  re- 
jquires  functions  two  orders  higher  than  four- 
jpoint  interpolation,  which  in  turn  was  two  orders 
[higher  than  the  functions  of  two-point  interpola- 
tion. Specifically,  six-point  interpolation  is  based 
on  a  sliding,  transhaping,  fourth -power  polyno- 
mial, two  orders  up  from  the  sliding  parabola. 
Each  arc  of  the  interpolation  locus  is  given  by  a 
fifth-power  polynomial. 

Six  mathematical  restraints  are  required.  In 
addition  to  the  four  used  for  four-point  interpola- 
tion, commonality  of  function  and  first  deriva- 
tive, commonality  of  the  second  derivative  with 
'the  second  derivative  of  the  linearly  transhaping, 
fourth-power  polynomial  is  required  at  the  ends 
of  the  interpolation  arcs.  The  junctures  of  the 
arcs  at  the  base  points  are  mathematically 
continuous  and  have  continuous  first  and  second 
[derivatives. 

Details  of  the  solution  of  the  coefficients  of  the 
six-point  method  are  not  repeated  here.  The 
solution  for  the  coefficients  in  equation  1  is 
conveniently  expressed  by  the  matrix-vector 
equation 
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In  the  previous  section  on  four-point  interpola- 
tion, it  was  pointed  out  that  the  method  could  be 
[extended  to  integration  of  experimental  data.  An 
advantage  of  the  six-point  system  is  that  it  can  be 
converted  to  differentiation,  leading  to  deter- 
mination of  smooth  rates  of  change  of  experi- 
mental data.  The  rate  of  change  of  values  in 
equation  1  is  given  directly  by  equation  2.  With 
coefficients  6  through  /  given  by  equation  43, 
rates  of  change  can  be  computed  by  direct 
substitution  in  equation  2.  An  advantage  of  this 
method  of  interpolative  differentiation  is  that  the 
first  derivatives  are  smoothly  continuous  be- 
cause the  second  derivatives  are  continuous. 


Figure  4  shows  six-point  interpolation  through 
and  the  calculated  first  derivative  of  the  same 
data  as  in  figure  3.  The  interpolative  values  and 
derivatives  were  computed  with  program  A-3 
(appendix  A). 

Six-point  interpolation  and  differentiation  can 
be  reduced  to  multiplying  coefficients  in  a  man- 
ner similar  to  interpolation  and  integration  with 
coefficients  as  presented  in  the  four-point 
method.  Interpolation  coefficients  for  each  one- 
hundredth  of  the  interval  of  the  interpolation  arc 
from  2=2  to  z=3  are  given  in  table  A-3  (appen- 
dix A).  These  coefficients  were  generated  by 
program  A-4.  Table  A-4  gives  coefficients  for  use 
in  differentiation.  Calculation  of  these  coeffi- 
cients is  included  in  program  A-4. 

It  should  be  noted  that  slopes  given  by  the 
differentiation  coefficients  of  table  A-4  are  in 
terms  of  the  scale  variate  z.  Given  that  the  scale 
of  original  data,  s,  is  related  to  z  by  s=kz-\-z0, 
then  ds=kdz.  A  slope  dy/dz  can  be  changed  to  a 
slope  of  original  data  by  simple  change  of 
variable,  yielding  dy/ds  =  (\/k)  dy/dz.  This  scale 
conversion  is  included  in  program  A-3,  but  must 
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FIGURE  4. — Six-point  interpolation  and  differentiation. 
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be  performed  as  a  separate  operation  following 
use  of  table  A-4. 

SMOOTHING  OF 
EXPERIMENTAL  DATA 

Smoothing  of  experimental  data  is  essentially 
the  application  of  the  method  of  least  squares  to 
effect  a  transformation  from  recorded  data  to  a 
set  of  base  points  for  interpolation.  Data  may  be 
recorded  at  nonuniform  intervals.  Transforma- 
tion produces  base  points  at  uniform  intervals. 
Data  may  include  mensuration  errors.  These 
errors  may  be  smoothed  by  averaging  during  the 
transformation.  Successful  averaging  requires 
that  base  points  be  fewer  in  number  than  data 
points  and  also  requires  a  configuration  of  data 
points  conducive  to  production  of  good  averages. 

The  location  and  spacing  of  the  base  points  is  a 
problem  in  design.  No  hard  rules  can  be  stated. 
The  heart  of  the  problem  is  the  amount  of 
averaging  that  is  desired.  With  a  goodly  amount 
of  data  in  a  compact  mass,  a  high  degree  of 
averaging  can  be  used.  The  base  points  would 
likely  be  one-fifth  to  one-tenth  of  the  number  of 
data  points.  With  limited  data,  strung  out  in 
range,  the  base  points  would  likely  be  about 
one-half  the  number  of  data  points,  resulting  in  a 
limited  degree  of  smoothing. 

In  the  following  examples,  a  high  degree  of 
smoothing  will  be  evident.  The  placement  of 
points  is  usually  dictated  by  the  particular 
problem.  A  base  point  would  usually  be  placed  at 
an  interface  or  boundary,  such  as  the  surface  of 
ground  or  water.  In  time  series  one  point  would 
usually  be  placed  at  time  zero.  Spacing  of  the 
remaining  base  points  through  the  data  range  is 
obviously  dictated  by  the  degree  of  averaging 
and  hence  the  number  of  points. 

Two-Point  Smoothing 

Two-point  smoothing  is  identical  to  evaluation 
of  connected  linear  segments  (Snyder  1967b). 
Equation  14  can  be  generalized  as 

y=C0YB0+Cl-YBl.  (44) 

In  equation  44,  the  meaning  is  limited  to  one 
interpolation  interval.  If  several  adjacent  inter- 
vals are  considered,  equations  45-47  result: 

y^COYBO+Cl-YBl  (45) 
y2~  C0YB1+C1YB2  (46) 

y3=  C0-YB2+C1YB3  (47) 


In  equation  45,  y  is  any  value  lying  on  the  line 
from  YBO  to  YBl.  In  equation  46,  y,  is  any  value 
lying  on  the  line  from  YBl  to  YB2,  and  so  on.  The 
interpolation  coefficients  CO  and  CI  are  known 
functions  of  scale  variate  z,  as  in  table  1. 

Assume  that  experimental  data  consist  of 
simultaneous  coordinate  values  of  x  and  y  for  a 
series  of  observations.  Then,  consider  a  series  of 
base  points  to  be  substituted  for  the  data  points. 
It  is  necessary  to  specify  the  location  of  the  base 
points  within  the  x-scale  of  the  data.  Normally, 
one  has  some  idea  of  the  number  of  base  points 
desired  to  span  the  data  in  the  x-scale.  Setting 
the  location  of  the  base  points  in  x-scale  automat- 
ically sets  a  2-scale  value  for  each  data  point  since 
the  base  points  are  always  spaced  at  unit  interval 
in  z. 

Referring  now  to  equations  45-47,  note  that 
they  are  linear  in  form.  Note  also  that  given 
experimental  values  of  y  and  given  a  computed 
e-variate  for  each  observation,  there  is  an  equiv- 
alent of  a  series  of  observational  equations  as 
used  in  multiple  regression.  CO  and  CI  take  the 
position  of  x0  and  x,  in  regression,  and  YBO  and 
YBl  take  the  position  of  regression  coefficients, 
conventionally  designed  b0  and  bx.  There  is  no 
regression  intercept  normally  designated  a.  It  is 
a  simple  matter,  with  conventional  least-squares 
methods  or  programs,  to  compute  YBO  and  YBl 
as  regression  coefficients.  A  transformation  from 
experimental  data  in  x-y  coordinates  to  base 
point  data,  in  z-YBl  coordinates,  is  thus  accom- 
plished. Note  especially  that  in  the  configuration 
of  equations  45-47,  YBl  and  YB2  must  take 
values  on  two  adjacent  intervals.  The  continuity 
of  the  smoothing  operation  is  thus  satisfied  with 
no  gaps  at  the  base  points. 

Table  B-l  (appendix  B)  shows  the  details  for 
organization  of  30  data  points  into  format  for 
regression  transformation.  The  data  points  and 
the  result  of  least-squares  transformation  in 
two-point  smoothing  are  plotted  in  figure  5. 
While  the  location  of  base  points  at  x=0.4,  0.6, 
0.8,  and  1.0  is  arbitrary,  the  values  of  these  base 
points  are  wholly  objective  and  the  result  of 
least-squares  fitting.  Note  specifically  in  figure  5 
the  overall  fit  of  the  linear  segments  to  the  rough 
crescent  shape  of  the  experimental  data.  Since 
the  two  internal  base  points  are  each  required  to 
be  common  to  two  segments,  the  fit  of  a  single 
segment  to  the  data  points  in  that  segment  is 
secondary. 
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FIGURE  5. — Smoothing  by  two-point  method. 
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FIGURE  6. — Smoothing  by  four-point  method. 


Four-Point  Smoothing 

Four-point  smoothing  is  a  simple  extension  of 
two-point  smoothing.  Equation  40,  modified  for 
adjoining  arcs,  is  extended  in  equations  48-51, 
analagous  to  the  two-point  equations,  45-47. 

y]=C0-YB0+Cl-YBl+C2-YB2+C3-YB3  (48) 
y  =  C0-YB2+C1-YB2+C2-YB3+C3-YB4  (49) 

2/3=  C0-YB2+C1-YB3 

+C2-YB4+C3-YB5  (50) 
y4  =  C0-YB3+C1-YB4 

+  C2-yB5+C3-YB6  (51) 


Values  for  the  coefficients  are  given  in  table 
A-l.  Remember,  in  equations  45-47  each  base 
point  is  linked  to  two  adjacent  intervals.  YBI  and 
YB2  in  those  equations  are  each  linked  by  CO  and 
CI.  Note  in  equations  48-51  that  an  internal  base 
point  such  as  YBS  is  linked  to  four  interpolation 
arcs.  The  least-squares  value  of  YB3  is  depend- 
ent on  experimental  data  in  four  intervals. 
Exterior  base  points  are  dependent  on  1,  2,  or  3 
intervals. 

Table  B-2  (appendix  B)  shows  organization  of 
experimental  data  in  format  for  four-point 
smoothing.  The  same  data  are  used  as  for  the 
two-point  smoothing  example,  except  the  x-scale 
is  shifted  to  the  left  to  put  the  first  base  point  at 
x=0.4.  The  z-scale  is  not  affected.  Program  B-l 
(appendix  B)  is  used  to  perform  the  least-squares 
transformation.  This  program  also  generates  the 
locus  of  interpolated  values  which  runs  through 
the  derived  base  points  and  generates  the  inter - 
grals  of  the  interpolation  arcs. 

Figure  6  shows  the  results  of  the  four-point 
smoothing.  An  interpolation  locus  of  smoothly 
joined  arcs  is  generated  that  runs  continuously 
through  the  data.  Comparison  of  figures  5  and  6 


shows  that  the  two-point  and  four-point  methods 
produce  base  points  not  greatly  different.  The 
main  difference  is  the  creation  of  arcs  by  the 
four-point  method,  smoothly  joined  with  angles 
eliminated.  Values  of  the  derived  base  points  can 
be  compared  in  table  2. 

The  integral  of  the  interpolation  arcs  is  also 
plotted  in  figure  6.  This  integral  can  be  thought 
of  as  the  smoothed  summation  of  the  experi- 
mental data. 

Six-Point  Smoothing 

Six-point  smoothing  is  a  simple  extension  of 
methodology  already  presented.  Equations  52-55 
are  analogous  to  equations  45-47  and  48-51. 

?/,  =  C0YB0+Cl-YB1  +  C2YB2 

+C3-YB3+C4-YB4+C5-YB5  (52) 
y2=C0-YBl+Cl-YB2+C2YB3 

+C3-YB4+C4-YB5+C5-YB6  (53) 
y  =C0-YB2+C1-YB3+C2YB4 

+C3  YB5+C4-YB6+C5  YB7  (54) 
y  =  CO  ■  YB3 + CI  ■  YB4 + C2 .  YB5 

+  C3YB6-I-C4YB7+C5-YB8  (55) 

In  equations  52-55,  note  that  an  interior  base 
point  will  depend  upon  values  of  experimental 


TABLE  2.  —  Comparison  of  derived  base  points 


Two-Point 

Four-Point 

Six-Point 

x  YBi 

X 

YBi 

X 

YBi 

0.4 

-4.14772 

0.4 

-1.72742 

.6 

-1.55935 

0.4  1.52418 

.6 

1.42889 

.8 

1.417502 

.6  3.75141 

.8 

3.59534 

1.0 

3.59252 

.8  3.81551 

1.0 

3.67674 

1.2 

3.68015 

1.2  2.57954 

1.2 

2.44932 

1.4 

2.45502 

1.4 

-  .42235 

1.6 

-2.29737 

1.8 

-20.17608 
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data  in  six  adjacent  intervals  by  linkage  through 
the  coefficients.  Exterior  base  points  may  de- 
pend on  five  or  fewer  intervals  of  data. 

Table  B-3  (appendix  B)  shows  organization  of 
experimental  data  in  format  for  six-point  smooth- 
ing. Program  B-2  is  used  to  perform  the  least- 
squares  transformation.  Following  transforma- 
tion of  the  experimental  data  to  a  smoothed  set  of 
base  points,  program  B-2  computes  a  locus  of 
interpolated  values  and  a  locus  of  differentials  of 
the  interpolation  arcs.  Figure  7  is  a  plot  of  the 
experimental  data  and  derived  interpolation  and 
differentiation  loci. 
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FIGURE  7. — Smoothing  by  six-point  method. 


Figure  7  contains  the  same  experimental  data 
as  figures  5  and  6.  It  can  be  seen  that  the 
interpolation  arcs  of  the  four-point  and  six-point 
methods  are  almost  identical.  Numerical  com- 
parison of  the  derived  base  values  is  shown  in 
table  2. 

Figure  7  also  shows  the  differentiation  of  the 
interpolation  arcs.  This  locus  of  values  can  be 
regarded  as  the  smoothed  rate  of  change  of  the 
experimental  data.  Note  the  smoothness  of  this 
locus  through  the  base  points  and  the  near 
constant  value  of  slope  from  about  x=1.14  to 
x=1.2. 

ATTRIBUTES  OF 
SLIDING  POLYNOMIALS 

Figures  5,  6,  and  7  demonstrate  that  various 
interpolation  methods  can  be  used  to  smooth 
experimental  data.  Smoothed  values,  whether 
read  from  two-point,  four-point,  or  six-point 


interpolation  loci,  are  not  extremely  different. 
Other  methods  would  also  produce  acceptable 
interpolation  loci.  Kimball's  (1974)  method  would 
require  a  means  to  separate  the  spectral  portions 
of  noise  and  information  in  the  data. 

The  identifiable  advantages  of  sliding  polyno- 
mials are  in  the  secondary  characteristics  of  their 
loci.  The  four-point  and  six-point  methods  are 
smooth  because  their  derivatives  are  continuous 
at  the  juncture  points  of  the  interpolation  arcs. 
The  four-point  method,  with  a  continuous  first 
derivative,  produces  an  extremely  smooth  locus 
of  integral  values,  since  integration  itself  is  a 
smoothing  process.  The  six-point  method,  with 
continuous  first  and  second  derivatives,  can  gen- 
erate a  smooth  rate  of  change  in  experimental 
data.  DuChateau  et  al.  (1972)  also  used  six  points 
in  evaluation.  In  their  method,  however,  para- 
bolic splines  are  required  to  be  continuous  only  in 
the  first  derivative.  Therefore,  these  splines 
cannot  be  expected  to  develop  the  smoothly 
continuous  rates  of  change  of  experimental  data. 
Such  continuous  smoothness  is  dependent  upon  a 
continuous  second  derivative. 

Sliding  polynomials  and  their  integrals  and 
derivatives  can  be  derived  from  experimental 
data  by  conventional  least-squares  techniques. 
The  interpolation  coefficients  provide  exact  error 
weights  for  each  experimental  data  point,  re- 
gardless of  the  interpolation  interval  in  which  the 
data  point  is  located.  The  structure  of  equations 
52-55,  for  example,  shows  that  the  evaluation 
procedure  is  similar  to  the  use  of  smoothing 
convolutes  as  in  Savitzky  and  Golay  (1964). 
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APPENDIX  A.  — COEFFICIENTS  AND  COMPUTER  PROGRAMS  FOR 

INTERPOLATION 

TABLE  A-l. —  Coefficients  for  four-point  interpolation 


1.00 

0.00000000 

1.OGOGG000 

•oooaoooo 

0. 00000000 

1.  01 
1.02 

.00490050 
.00960400 

. 99975450 
.99901200 

. 00519850 
.01078800 

-. 00004950 
-.00019600 

1.03 
1.04 

.01411350 
.01843200 

.99779050 
.99609600 

.01675950 
.02310400 

-.00043650 
-.00076800 

1.05 
1.06 

.02256250 
.02650800 

.99393750 
.99132400 

. 02981250 
.03687600 

-.00118750 
-.00169200 

1.07 
1.08 

.03027150 
.03385600 

. 98826450 
.984768U0 

.04^+28550 
.05203200 

-.00227850 
-.00294400 

1.09 

.03726450 

. 98084350 

.06010650 

-.  00  368550 

1.  10 

.04050000 

. 97650000 

. 0685u  00  0 

-. 00450C00 

1.11 
1.12 

.04356550 
.04646400 

.97174650 
.96659200 

. 07720350 
.08620800 

-.00538450 
-.00633600 

1.13 
1.  14 

.04919850 
.05177200 

•961C4550 
.95511600 

.09550450 
.10508400 

-.00735150 
-.  00  842800 

1.15 
1.16 

.05418750 
.05644800 

. 94881250 
.94214400 

.11493750 
.125G560Q 

-.00956250 
-.01075200 

1.17 
1.18 

.05855650 
.06051600 

.93511950 
.92774800 

.13543050 
.  14605200 

-.01199350 
-.01328400 

1. 19 

.06232950 

.92003850 

.15691150 

-.01462050 

1.20 

. 06400000 

.91200000 

•16800000 

-.01600000 

1.21 
1.22 

.06553050 
.06692400 

.90364150 
. 89497200 

. 17930850 
.  19082800 

-.01741950 
-.01887600 

1.23 
1.24 

.06818350 
.06931200 

. 88600050 
. 87673600 

.20254950 
.21446400 

-.02036650 
-.02188800 

1.25 
1.26 

.07031250 
.07118800 

. 86718750 
.85736400 

.22656250 
. 238836C0 

-.02343750 
-.02501200 

1.27 
1.28 

.07194150 
.07257600 

. 84727450 
. 83692800 

. 25127550 
.26387200 

-.02660850 
-.02822400 

1.29 

.07309450 

.82633350 

.27661650 

-.02985550 

1.30 

.07350000 

.81550000 

.2895000  0 

-.03150000 

1.31 
1.32 

.07379550 
.07398400 

•  804<*3650 
.793152C0 

.30251350 
. 31564800 

-.03315450 
-.03481600 

1.33 
1.34 

.07406850 
.07405200 

.78165550 
.76995600 

. 32889450 
.3422440  0 

-. 03648150 
-.03814800 

TABLE  A-l. — Coefficients  for  four-point  interpolation — Continued 


1.35 
1.  36 

— 

.07393750 
. 07372800 

.758C6250 
.74598400 

. 35568750 
.36921600 

.03981250 
. 04147200 

1.  37 
1.38 

.07342650 
.07303600 

.73372950 
.72130800 

.38282050 
. 39649200 

.04312350 
. 04476400 

1.39 

• 

.07255950 

.70872850 

.41022150 

.04639050 

1.40 

.07200000 

.69600000 

.42400000 

.04800000 

1.41 
1.  43 

.07136050 
.07064400 

.68313150 
.67013200 

.43781850 
. 45166800 

.04958950 
. 05115600 

1.  43 
1.44 

.G6985350 
.06699200 

.65701050 
.64377600 

.46553950 
.479*2400 

. 05269650 
. 05420800 

1.45 
1.46 
1.47 
1.  48 

.06806250 
.06706800 
.06601150 
. 06489600 

.63043753 
. 61700400 
.60348450 
.56988600 

.49331250 
.50  719600 
.52136550 
.53491200 

-.05568750 
-.05713200 
-.05853850 
-. 35990400 

1.49 

.06372450 

.57622350 

.54872650 

-.06122550 

1.50 

.  0  6250  300 

•  562 50 OC  J 

.5625030  0 

-.06250000 

1.51 
1.52 

.06122550 
.05990400 

.54872650 
.53491200 

.57622350 
.58988800 

-.06372450 
-.06489600 

1.53 
1.  54 

.05853850 
.35713200 

.52106550 
.50719600 

.60348450 
.617jG*G0 

-. 06601150 
-.06706800 

1.55 
1.56 

.05568750 
.05420800 

.49331250 
.4794240C 

.63043750 
.64377600 

-. 06806250 
-.36899200 

1.57 
1.58 

.05269650 
.05115600 

.46553950 
.45166800 

.65701050 
.67013200 

-.06985350 
-.07064400 

1.59 

.04958950 

.43781850 

.68313150 

-.07136050 

1.60 

.04800000 

.42*00000 

•6960uu00 

-.07200000 

1.61 
1.62 

.04639350 
.04476400 

.41022150 
. 3964920  0 

.70872850 
.72133800 

-.07255950 
-.073036C0 

1.63 
1.64 

.04312350 
.04147200 

.38282050 
.36921600 

.73372950 
.74598400 

-.07342650 
-.07372800 

1.65 
1.66 

-.03981250 
-.  0  3814800 

.35568750 
. 34224400 

.75806250 
.76995600 

-.07393750 
-.07405200 

1.67 
1.68 

-.03648150 
-.03481600 

.32889450 
. 31564800 

.78165550 
.7931520  0 

-.07406850 
-.07398400 

1.69 

-.03315450 

.30251350 

.804*3650 

-.07379550 

TABLE  A-l.  —  Coefficients  for  four-point  interpolation— Continued 


1.70 

•03150000 

. 28950000 

•  81550000 

-.07350000 

1.71 
1.72 

.02985550 
.028221*00 

.27661650 
. 26387200 

.82633350 
. 83692800 

-.073091*50 
-.07257600 

1.73 
1.  7k 

.02660850 
.02501200 

.25127550 
. 23883600 

.81*7271*50 
.  857361*00 

-.07191*150 
-.07118800 

1.75 
1.76 

.0231*3750 
.02188800 

. 22656250 
.  21^6<*Q0 

.86718750 
. 87673600 

-.07031250 
-. 06931200 

1.77 
1.78 

.02036650 
.01887600 

. 20251*950 
.  19082800 

.88600050 
. 89U97200 

-.06818350 
-.066921*00 

1.79 

.0171*1950 

. 17930850 

.90361*150 

-.06553050 

1.80 

.01600000 

.16800000 

.912G0000 

-.  061*00000 

1.81 
1.  82 

.011*62050 
.  0  13281*00 

.15691150 
.  11*605200 

.92003850 
.  92771*800 

-.06232950 
-. 06051600 

1.  83 
1.81+ 

.01199350 
.01075200 

. 1351+3050 
.12505600 

. 93511950 
.  9^*211*1*00 

-.05855650 
-.0561*1*800 

1.85 
1.86 

.00956250 
•  008<+2800 

.111*93750 
. 105081+00 

.91*881250 
.95511600 

-.  051*18750 
-.05177200 

1.87 
1.88 

.00735150 
.00633600 

.095501*50 
.08620800 

.96101+550 
. 96659200 

-.0U919850 
-.  0i*6i*6i*00 

1.89 

.005381*50 

.07720350 

.97171+650 

-.01*356550 

1.90 

.001*50000 

.06850000 

.97650000 

-.  01*050000 

1.91 
1.92 

.00368550 
.0  0291*1*00 

.06010650 
.05203200 

.98081*350 
.981+76800 

-.037261*50 
-.03385600 

1.93 
1.9V 

.00  227350 
.  0  0 1 69  2  0  0 

.0*+ 1*28550 
.03687600 

.988261*50 
.991321+00 

-.03027150 
-.02650800 

1.95 
1.96 
1.97 
1.98 

.00118750 
.00076300 
.  0001*3650 
•U00196UO 

.02981250 
.  023101*00 
.01675950 
.01078800 

.99393750 
.99609600 
. 99779050 
.99901200 

-.02256250 
-.  0181+3200 
-.  011+11350 
-.  009601+00 

1.99 

.00  0  01*950 

.00519850 

.99975150 

-.  001*90050 

2.  30 

C .OuQOOOOG 

-.00000000 

1. 00000000 

0. 00000000 

TABLE  A-2. —  Coefficients  for  four-point  integration 


14? 

CO 

CI 

C2 

C3 

,.,  ?•     „   V  fl  ; 

1.00 

- 

.00000000 

-. Q00G0000 

O.Oi-OuQOQO 

.00000000 

4      ft  4 

1  •  U  1 

1.02 

•  U  U  U  Uc  *4D  f 

.00009735 

nno QQQ1 7 
.01999339 

.  u  u  u  u  c  ?b  o 
.  00  010527 

_    n  n  ft  ft  n  n  4  7 
-.00000131 

1.0  3 
1.04 

.00021610 
.00037899 

.02997780 
.0  3994763 

.00024270 
.00044171 

-. 00000440 
-.00001035 

1.05 
1.  06 

.00058411 
.00082962 

.04989818 
. 05982486 

.00070599 
. UG 103914 

-. 00002005 
-.00003438 

1.07 
1.08 

.00111367 
.00143445 

.06972317 
. 07958869 

.00144466 
.00192597 

-.00005417 
-.00008021 

1.09 

.00179020 

.08941710 

.  00  248640 

-. 00011330 

1.10 

.00217917 

.09920417 

.00312917 

-.00015417 

1.11 
1. 12 

.00259963 
.  00  304992 

. 10894574 
. 11863776 

.  00  385743 
.00467424 

-.00020353 
-.00026208 

1.13 
1.14 

.00352837 
.00403335 

.12827627 
.13785739 

.00558256 
.00658527 

-.00033047 
-.00040931 

1.15 
1. 16 

.00456328 
.0U511659 

. 14737734 
. 15683243 

•C0768516 
. 00  888491 

-.00049922 
-.00060075 

1.17 
1.18 

.00569173 
.00628722 

.  1662190'* 
.17553366 

.01018713 
.01159434 

-.00071443 
-.00084078 

1.19 

.00690157 

. 18477287 

.01310896 

-.00098027 

1.20 

.00753333 

.19393333 

.01473333 

-.00113333 

1.21 
1.22 

. 00818110 
.00884349 

.20301180 
.21200513 

.01646970 
.01832021 

-.00130040 
-.00148185 

1.  23 
1.24 

. 00951913 
.01020672 

. 22091024 
.22972416 

. 02028693 
.02237184 

-. 00167803 
-.00188928 

1.25 
1.26 

. 01090495 
.01161255 

. 23844401 
.24706699 

.02457682 
.  02690  367 

-. 00211589 
-.00235811 

1.27 
1.28 

.01232830 
.01305099 

. 25559040 
.26*01163 

.02935410 
.03192971 

-.00261620 
-.00289035 

1.29 

.01377943 

. 27232814 

.03463203 

-.00318073 

1.30 

.01451250 

.28053750 

.03746250 

-.00348750 

1.31 
1.32 

.01524907 
.01598805 

.28863737 
. 296625*9 

. 04042246 
. 04351317 

-. 00381077 
-.00415061 

1.33 
1.34 

.0167281*0 
.0171*6909 

.30449970 
. 31225793 

. 04673580 
.05009141 

-.00450710 
-.00488025 

TABLE  A-2.—  Coefficients  for  four-point  integration  —  Continued 


1.  35 
1.36 

- 

.01820911 
.0 1894752 

. 31989818 
.32741856 

.05 35 8099 
.05720544 

-.005 27 005 

•    WW        U>  9    W  W  ^ 

-.00567648 

1.37 
1.  38 

- 

.01968337 
.02041575 

.33481727 
. 34209259 

.06096556 
.66486207 

-• 00609947 

•     WW    W  W    J    J  ' 

-.00653891 

1.  39 

.02114380 

. 34924290 

•  06  889560 

-. 00699470 

1.  40 

- 

.02186667 

. 35626667 

. 07306667 

-.00746667 

1.41 
1.42 

.02258353 
.02329362 

.36316244 
. 36992  886 

. 07737573 
.08182314 

-. 00795463 
-. 00  845838 

1.43 
1.44 

.02399617 
.02469045 

. 37656467 
. 38306869 

.08640916 
.09113397 

-.00897767 
-.00951221 

1.45 

.02537578 

. 38943984 

.09599766 

-.01006172 

1.46 

.02605149 

. 39567713 

. 10100021 

-. 01062585 

1.47 

.02671693 

.40177964 

.1G61«*153 

-. 01120423 

1.48 

.02737152 

. 40774656 

.11142144 

-.01179648 

1.49 

.02801467 

.41357717 

.11683966 

-. 01240217 

1.50 

.02864583 

• 

. 41927083 

.12239583 

-.01302083 

1.51 

.  Q2926<t50 

.42482700 

. 12808950 

-. 01365200 

1.52 

.62987019 

.439  24523 

.13392011 

-.01429515 

1.53 

.03046243 

.43552514 

.1398870  3 

-.01494973 

1.54 

.03104382 

.44066646 

.14598954 

-.01561518 

1.55 
1.56 

.0  3160  495 
.0  3215445 

.44566901 
.45053269 

. 15222682 
.15859797 

.01629089 
. 01697621 

1.57 
1.58 

.03268900 
.0  3320  329 

.45525750 
.45984353 

. 16510200 
.17173781 

. 01767050 
.01837305 

1.59 

.03371203 

.46429094 

.17850423 

. 01908313 

1.60 

.0  3420000 

.46860000 

.18540000 

.01980000 

1.61 
1.62 

.03467197 
.03512775 

.47277107 
.47680459 

. 19242376 
.19957407 

. 02052287 
. 02125091 

1.63 
1.64 

.03556720 
.0  3599019 

.48070110 
.48446123 

.20684940 
.21424811 

. 02198330 
.02271915 

1.65 
1.66 

.0  3639661 
.03678642 

.48808568 
. 49157526 

.22176849 
.22940874 

.02345755 
.02419758 

1.67 
1.68 

.03715957 
.03751605 

.49493087 
.49815349 

.23716696 
.24504117 

. 02493827 
. 02567861 

1.69 

.03785590 

.5C124420 

.25302930 

.02641760 
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TABLE  A-2.  —  Coefficients  for  four-point  integration  —  Continued 


1.70 

.0  3817917 

•  5042041 7 

Of    A  A    O  f\  A 

•  26112917 

ft  O  ^9  A   r*  i    A  —J 

- • 02715417 

1.71 
1.  72 

- 

.03848593 
•  0  3877632 

.50703464 
•  50973696 

.26933853 
•27765504 

-.02788723 

n  oi  a  f  a  ^  f  a 

- • 0  2861568 

1.73 
1.74 

.03905047 
.0  3930855 

.51231257 
.51476299 

.28607626 
.29459967 

-.02933837 
-.03005411 

1. 75 
1.76 

- 

•  0  3955 u  7 8 
.03977739 

^m  A  *9  t~    ft        ft  j 

•  5 171 8984 
.51929^83 

•  3  0  322266 
.31194251 

ft  — *  ft  4r    A  — •  A 

-. 03076172 
-.03145995 

1.77 
1.  78 

. 0  3998863 
.04018482 

. 52137974 
. 52334646 

• 32  075643 
. 32966154 

ft  ^  O  A   t    T  V  *9 

- .  0  3214753 
-. 03282318 

1.79 

.04036627 

.52519697 

.33865486 

-.03348557 

1.80 

.04053333 

. 52693333 

. 34773333 

-. 03413333 

1.  81 
1.82 

.  U4G68640 
.04082589 

.5285577U 
.53007233 

" .3568938  0 
.36613301 

-.03476510 
-.03537945 

1.83 
1.84 

.04095223 
.04106592 

.53147954 
. 53278176 

. 37544763 
.38483424 

-.03597493 
-.03655008 

1.85 
1.86 

.04116745 
.04125735 

.53398151 
.53508139 

.39428932 
.40380927 

-.03710339 
-. 03763331 

1.87 
1.88 
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EXAMPLES  USING  TABLES  A-l  AND  A-2 


EXAMPLES  USING  TABLES  A-3  AND  A-4 


Example  1 .  — Interpolate  a  value  of  Y  at  x=0.07 

}iven  the  table  of  base  points 


X 

BYi 

-0.05 

1.00 

.00 

1.15 

.05 

1.25 

.10 

1.25 

.15 

1.10 

.20 

.75 

.25 

.15 

Example  3.  — Given  the  same  base  points  as  in 
example  1  interpolate  at  x=0.07. 

From  table  A-3: 


(C0=0 

0144)xl.00= 

0.0144 

(Cl=- 

.11361x1.15= 

-.1306 

(C2= 

.7264)xl.25= 

.9080 

(C3= 

.43841x1.25= 

.5480 

(C4=- 

.0736)xl.l0  = 

-.0810 

(C5= 

.0080)  x  .75= 

.0060 

X  0.07  = 

1.2648 

z=l  +  (0.07-0.05)/(0.10-0.05)=1.4 


?rom  table  A-l: 

(C0=  -0.072)xl.l5=  -0.0828 

(Cl=    .6961x1.25=  .8700 

(C2=    .424)xl.25=  .5300 

(C3=  -.0481x1.10=  -.0518 

V  o.o7  =  l-2644 


Example  2.  — Compute  the  integral  from  x  =  0.00 
ox  =  0.07 


Example  4.  — Interpolate  for  the  differential  of  y 
at  x  =  0.07 

From  table  A-4: 

(C0=0. 018001x1. 00=  0.0180 

(C\=  .097331x1.15=  .1119 

(T2=l. 322671x1. 25=  1.6533 

(T3=  1.464001x1.25=  1.8300 

(C4=  .259331x1.10=  -.2853 

~(C5=  .03867)x  .75=  .0290 

.0143 


kly/dx)  „  ()7  =(0.0143/0.05)  =  0.286 


from  0.00  to  0.05 
rrom  table  A-2: 

(C0=-0.04167)xl.00  = 
(Cl=  .541671x1.15= 
(C2=  .541671x1.25= 
(C3=  -.041671x1.25: 


From  0.05  to  0.07 
?rom  table  A-2: 

(C0=-0. 02187)xl. 15: 
(Cl=  .356271x1.25: 
(C2=  .073071x1.25: 
(C3=  -.007471x1.10: 


NOTES  FOR  PROGRAMS  A-l,  A-2,  A-3,  AND  A-4 


-0.0417 
.6229 
.6771 
-.0521 
1.2062 


-0.0252 
.4453 
.0913 

-.0082 

.5032 


/ 

•/a  i 


0.07 


0.00 


ydx  =(1.2062  +  0.50321x0.05  =  0.08547 


TAB(I) 


YINT(I) 
YAINT(I) 
YI'INT(I) 
INTFORl 

INTSIX( 
CHI) 

ICOEF(I) 

YAO(I) 
XKI) 


Input  of  tabular  values,  or  base 

points,  on  which  interpolation  is 

based. 
The  interpolated  values. 
Interpolated  integral. 
Interpolated  derivative. 
)  Subroutine  for  four-point 

interpolation. 
)    Subroutine  for  six-point 

interpolation. 
Coefficients  of  equations  37  and  43. 
The  matrices  of  equations  37  and 

43. 

Coefficients  of  equation  42. 
The  index  z  of  tables  A-l,  A-2,  A-3, 
and  A-4. 
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PROGRAM  A-1.  —  Four-point  interpolation  and  integration 


DIMENSION   YINT(20G ) , YAINT (200 ) 
COMMON/XX/TAEM1Q0 ) 
READ (5, 10 00) (TA3 (I) ,1=3,22) 
10  0  0   FORMAT  ( 20  F4 • 2) 
TAP  ( 1 )  =  0 . 1 
TAR (2) =0. 1 

X  =  0.3  ___™ 
MS=-1 

no  iooi  1=1,151 

X=X+0.31 

CALL   INTFOR (X,Y ,YA ,MS) 

YINT(I) =Y 

Y  AINT ( I )  =  Y  A 
1001  CONTINUE 

WRITE (6, 10 02) (I, YI NT ( I) ,1=1,151) 
1  00  2  CCRMAT  (  *   INTERPOLATED  FUNCTION*/ (8  (15,  F10  .<♦))  ) 

WRITE (6, 100 3)    ( I, Y AINT (I) ,1=1,151) 
1003  FORMAT (*   INTERPOLATED  INTEGRAL*/ (8 (15, F 10 . k) ) ) 

STOP 

END 


SUBROUTINE   INTFOR ( X, Y ,YA , MS) 
DIMENSION   CCf)  ,  ICOEF  (4,4) 
C0MM0N/XX/TAD( 100 ) 

DATA   ICOEF/^, -8, 5, -1,-6, 19, -14, 3, 6, -16, 13,-3,-2,5,-^,1/ 
AX=( (X-TAB(l) ) /TAB (2 ) ) 
M  =  I NT ( A  X) 

F  X=  AX-FLO  AT  (M)  +1.0  

IF(M.EQ.MS)    GO   TO  102 

DO  100  1=1,4   

C(i)=0.Q 

DO   100  J=l,4 

100  C(I)=C(I)  +  ICOEF (I , J) *  TAB (M  +  J+l) 

CA  =  C(i)+C(2)/2.0  +  C(3)/3.0*-C(4)/4.0 
MS  =  ^ 
102  Y=C(4) 

Y A=C(4) /4. 0 
DO   101  K=l,3 
Y=Y*FX+C( 4-K) 

101  YA=YA*FX-»-C  (4-K)  /  <4-K)  _____ 
Y=Y/2. 3 

Y A= (YA*FX-CA) *TAB( 2) /2.0 

RETURN 
END 
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PROGRAM  A-2.— 


Generation  of  four-point  coefficients 


 Time n s i on"  y in?  ( 1  o  1 ,  4 ) ,  y  a  i n t  ( 1  o  i ,  4 > ,  1  co  e f  ( 4 ,  4 ) , x  1  ( 1 0 1  > ,  y a 0  ( 4 ) 

D  AT A  130  E  F  /  4 j^_8_,  5, ^l,-6jl 9, -14, 3, 6, -16, 13,-3,-2,5,-4,1/ 
Y AO ( 1) = -17. 0/24.3 

Y  AO  (2_)  =5.  0/24^0 

Y  AO (3 ) = - 1 9.0/24.0 

Y AO  (_4)_  =  7._0/ 24_._ 0__  _ 

"90 iOGl  J =1,101 

X  =  FLOAT  (Ji/lOO-JtG   

XI  ( J)=X 

90  10  Q 1  1  =  1,4   

1 1NK  J,  I )  =  ICOEF  ( 4  , 1 ) 
DO   10  0  2  <=J_l1„ 
1 0  0  2  Y  IN "f  <  J  >  I)  =  YI N T  ( j ~,  I )  * X  + 1 C 0 1 E F  ( 4 -  K  ,  I ) 

Y  I N  T  (  J, I)  =  Y I  NT  (J jj  I  )  /  2  .  0 

Y  AT. NT  ( J  ,  I  )  =  I COE  F  (4  , 1 )  /4 .  0 
3O__10  0  5  _K  =  1,  3 

1003  YAINT(w,I)=(YAINT (  j ,  I )  *X  )  +  <  FL  0  A T  <  I COEF <4-K,I) )/ (4-  K) ) 

Y  A I NT ( J , I  YAINT(J,I)/2.0)*X+YAO(I) 
1 0  0 1  C  ON  T IN U  E  ~ 

1  00  0  CONTINUE   _   

DO  10  0  5  J= 1,101 

IF(  J_.E3_.l  .OR.  J.EQ.  5_l_)_GO  TO  1006 
"GO  TO  100  7 

1006  WRITE  (6, 1C0 8Ji  _  

100  8  FORMAT (*i*> 

1  0  0  7  I  F  ( MOO  (J ,  51 .  EQ.  1 1  GO  TO  _100  9 
GO~  TO  101 C 

10  0  9  W_R_I  TE  (  6_i  1  0  11)     _ 

foil  FORMAT  1*0*-) 

1010   W  RITE ( o , 1 C  12 )  (XI ( J  ) ,  <  YINT ( J ,1 )  ,1=1,4)) 
1 0 1 2  F ORMAT (*   *  ,  (  35X  ,FE> .  2  ,4F12  .  8)  ) 

1  00  5  CONTINUE  _ 

D O  10  2 0  J  =  l,101 

I F ( J • E3 . 1 .OR .J.EQ. 51)   G O  TO  1 01 3 
GO  TO  1014 

1  0 1 J  W^H^_( 3  ll  P_ML  __  „  

1 0 1 4  1 F ( MOO ( J, 5 ) .tO • 1)    GO   TO  1016 

GO  TO  1017  

1016  WRITE  (6  ,1011)" 

1  0  17  H  RIT  E _{o  »!_?  1 2  >  <  *I  <  jJ  >  *  J  Y  ^  T  ( J,  I )  ,  I  =1 ,  4 )  ) 
10  20  CONTINUE 

STOP   _  _    

END  ~ 
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PROGRAM  A-3.  —  Six-point  interpolation  and  differentiation 


DIMENSION   YINT(20C) , YPINTC200) 
COMMON/XX/TAPttOO) 
REA0C5  ,1Q00MTAB<I1  ,1  =  3,22) 
1  00  0  FORMAT { 20 F*. 2) 
TAB(1)=0#1 

TA9(2)=0.1   

X  =  0.  3 

MS  =  -i   

DO  1001  1=1,151 
X  =  X«-0 . 3 1 

CALL   INTSIX(X,Y ,YP,MS> 
YINT(I) =Y 
YPINT(I)= YP 

1001  CONTINUE 

WRITE(6»1302)  tl » YINT  (II  ,  1= 1  ,151) 

1002  FORMAT  i*   INTERPOLATED  FUNCTION*/ (R ( 15, F 10 . k) )  ) 
WRITE (6, 100 3)  (I,YPINT  (I)  ,1=1,151) 

1  00  3  FORMAT ( *   INTERPOLATED  DE RI V AT  I VE*/ ( 8 (1 5 , F10 . k ) }> 
STOP 
END 


SUBROUTINE   INTSIX < X , Y , YP ,MS) 
DIMENSION     CT6) tIC0EFC6t 61 
CCMMON/XX/TAR(130 ) 

DATA  ICOEF/it32,-c'lfl,  765,  -  313,  6  3,-5  ,-20<*Q,  *i*36,- 375  k,  1551,-31^,25, 
1*080, -8 75?, 7 klk, -3078, 62 6, -50, -408T, 87 12, -73 5 6, 3058, -62*, 50, 2  0*0, 
2-^3^6,3661,  -1521,  311  ,-25  ,-<+0  8  ,  8  68,-730  ,  30  3,-6  2,5/ 

AX= ( (X-TAB(l) ) /TAP  (2) ) 

M=INT( AX) 

FX=AX-FLOAT(M) +2. 0 
IFCM.EQ.MS)   CO  TO  102 
DO   100  1=1,6 
C (I)=C  .  C 
0  0   1C0  J=l,6 
10  0   C  (I)=C( I) + ICOEF (I , J) *TAB(M+ J) 
MS=M 

10?  Y=CC6) *FX+C(5) 
Y  P  =  5  .  <3  *  C  (  6 ) 

no  101  K-lyk 

Y=Y*FX+C( 5-K) 
101  YP=YD*FX+ (5-K? *C ( 6-K) 
Y=Y/2U . 0 

YP=YP/2<*.  0  /TA8(2) 

RETURN 

END 
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PROGRAM  A-4.  —  Generation  of  six-point  coefficients 


DIMENSION   YINT(101,6)  f  YPINTdOl  ,6)  »IC0EF<6,6)  fXI  (101) 

nATA_.I0.0EF/432,- 91  8_,  765,  -  31  3,  6  3  ,  -5  ,-20  *»0  ,  k*36  3754  ,15  51 , -31V,  25*. 
l<+08C,-e75  2,7<+li+,-3078,526,-50,-i+C8C,87i2,-735  6,3C58,-6  2iM  50,20<»0, 
2-^346,366Jl,-i521,3.1l4--25,-<t0  8  ,  86  8, -7  30  ,  30  3  , -6  2,  5/ 

00   1000  J=l,101 
„X  =  EUQA.T  (J )  V1HQ_  Jltl *  99 

XI ( J) =x 
  0-0  1001  1  =  1,6  .  

YINT(J,I) =ICOEF (6,1) 

  OQ.  10  0  2  K=.l,5.   

100  2  YINT( J, I) =YINT( J , I ) *X +ICOE F (6 - K , I) 

YIMT  (J  ».I)  _=YIMT  {  J.,Ii/2*t.!J 

Y  PINT  (J,  I)  =IC0EF(6  ,1)  *5.  0 

0  0.  10  Q  3  K_=l,4 
LOO  3  Y  PI  NT  <(  J,  I)  =Y^INT(  J,  I  )  *X+ ICOEF ( 6- K , I ) * ( 5-K ) 

YPINT  tJ  ,I)-YPIN.T.(  J  ,1)  /2t*.0 
L001  00NTINUc 
LOGO  CONTINUE.   ..   

TO  1005  J=l,101 

I-F  (U*£3.1.nR.  J.  EQ,  51)    GO  TO  1006 

0  0  to   1Q0  7 

L0  06  WRITE (6 ,1008) 
1008  F0RMflT(*l*) 

L007  I  F( MOO (J, 5). ECU  IX  GO   TO  1C09 

GO  TO  101C 
10Q9  WRITE(6,10il) 
1011  FORMAT  C *0 *) 

1010  WRITE  (5,1012)  (XIU)  ,  (YINT  ( J,I)  ,1=1,6)) 
1312  FORMAT ( *   *, (35X,F5.2,6F12.8)) 

100  5  CONTINUE 

DO   1020  J=l,101 

1  F. ( J  .£.1. 1  .  OR .  J .  E.Q. .  5.1 )   GO  TO  1013 
GO  TO  101^ 

1013  WRITE(6,lC0ft)   _ 

13H+  IF(M0D(J,5)  .FQ.  1)    GO   TO  1016 

GO  to   131 7 
1016  WRITE (5,1011) 

101?  WRITE  <6_,1 012>  (XI  (J)  ,  (YPINT  (J,  I )  ,  1=1 ,6)  > 
1020  CONTINUE 

GTOP 

rNO 
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APPENDIX  B.  — DATA  ARRANGEMENT  AND  COMPUTER  PROGRAMS  FOR 

SMOOTHING 


TABLE  B-l.  -  Data  arrangement  for  two-point  smoothing 
red  Base-point         Scale  Values  of  coefficients 1 


data  points  positions         variate  for  base  points 


y 

in  x 

z 

YBl 

YB2 

YB3 

YB4 

V 

0  A 
U.4 

U.41Z 

1  so 

1 .0\J 

U.Uo 

0.94 

0.06 

1.50 

A  1  £ 
.410 

1  fiQ 
1  .Uf 

no 
.U5 

.92 

.08 

1.69 

a 

.405 

9  09 

OO 

.zy 

.71 

.29 

2.02 

a  i^n 

.40U 

9  90 

.75 

.25 

2.20 

AP.A 

9  Sfl 

.66 

.68 

.32 

2.50 

.ouu 

9  50 

.ou 

.50 

.50 

2.50 

^oo 
.ouu 

2.70 

.ou 

.50 

.50 

2.70 

.514 

2.89 

.57 

.42 

.57 

2.89 

.566 

3.25 

.83 

.17 

.83 

3.25 

R1A 
.0/4 

O.JO 

.O  1 

1  O 

.13 

.87 

3.58 

.ouu 

3  fi9 

1  no 
1  .UU 

1.00 

3.62 

.0^4 

o .  oo 

1  o 

.88 

0.12 

3.66 

.  o^o 

3.85 

1  9. 

.87 

.13 

3.85 

.oou 

^  s^ 

.30 

.70 

.30 

3.55 

.oou 

4.U4 

.30 

.70 

.30 

4.04 

.000 

o.oo 

A  9 

.43 

.57 

.43 

3.88 

.oyo 

3  fiO 

.4y 

.51 

.49 

3.60 

.  I  10 

^  Q8 

CO 

.oo 

.42 

.58 

3.98 

.750 

3.70 

.75 

9S 
.£0 

7S 

3.70 

.774 

3.87 

.87 

.13 

.87 

3.87 

.806 

3.58 

.03 

.97 

.03 

3.58 

.844 

3.69 

.22 

.78 

.22 

3.69 

.860 

3.29 

.30 

.70 

.30 

3.29 

.866 

3.56 

.33 

.67 

.33 

3.56 

.886 

3.27 

.43 

.57 

.43 

3.27 

.910 

3.10 

.55 

.45 

.55 

3.10 

.914 

3.23 

.57 

.43 

.57 

3.23 

.956 

2.74 

.78 

.22 

.78 

2.74 

.960 

2.98 

.80 

.20 

.80 

2.98 

.986 

2.55 

.93 

.07 

.93 

2.55 

1.0 


^1        X2        X3        X4  Y 
Conventional  regression  symbols 

Equations  12  and  13. 
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TABLE  B-2.  —Data  arrangement  for  four-point  smoothing 


Base-point  Scale  Values  of  coefficients1 
positions         variate   for  base  points 


in  x 

z 

YBl 

YB2 

YB3 

YB4 

YB5 

YB6 

Y 

0.4 

.6 

1.06 

-0.0265 

0.9913 

0.0369 

-0.0017 

1.50 

1.08 

-.0339 

.9848 

.0520 

-.0029 

1.69 

1.29 

-.0731 

.8263 

.2766 

-.0299 

2.02 

1.25 

-.0703 

.8672 

.2266 

-.0234 

2.20 

1.32 

-.0740 

.7932 

.3156 

-.0348 

2.50 

1.50 

-.0625 

.5625 

.5625 

-.0625 

2.50 

1.50 

-.0625 

.5625 

.5625 

-.0625 

2.70 

1.57 

-.0527 

.4655 

.6570 

-.0699 

2.89 

1.83 

-.0120 

.1354 

.9351 

-.0586 

3.25 

1.87 

-.0074 

.0955 

.9610 

-.0492 

3.58 

2.00 

.0000 

.0000 

1.0000 

.0000 

0.0000 

3.62 

1.12 

-.0465 

.9666 

.0862 

-.0063 

3.66 

1.13 

-.0492 

.9610 

.0955 

-.0074 

3.85 

1.30 

-.0735 

.8155 

.2895 

-.0315 

3.55 

1  .ou 

.1/  1  oo 

4.04 

1.43 

-.0699 

.6570 

.4655 

-.0527 

3.88 

1.49 

-.0637 

.5762 

.5487 

-.0612 

3.60 

1.58 

-.0512 

.4517 

.6701 

-.0706 

3.98 

1.75 

-.0234 

.2266 

.8672 

-.0703 

3.70 

1.0 

1.87 

-.0074 

.0955 

.9610 

-.0492 

3.87 

1.03 

.0000 

-.0141 

.9978 

.0168 

-.0004 

3.58 

1.22 

-.0669 

.8950 

.1908 

-.0189 

3.69 

1.30 

-.0735 

.8155 

.2895 

-.0315 

3.29 

1.33 

-.0741 

.7417 

.3289 

-.0365 

3.56 

1.43 

-.0699 

.6570 

.4655 

-.0527 

3.27 

1.55 

-.0557 

.4933 

.6304 

-.0681 

3.10 

1.57 

-.0527 

.4655 

.6570 

-.0699 

3.23 

1.78 

-.0189 

.1908 

.8950 

-.0669 

2.74 

1.80 

-.0160 

.1680 

.9126 

-.0640 

2.98 

1.93 

-.0023 

.0443 

.9883 

-.0303 

2.55 

1.2 
1.4 


X\       X<i       X3       X4       X5       Xq  y 
Conventional  regression  symbols  


1 

From  table  A-l. 


TABLE  B-3.  -Data  arrangement  for  six-point  smoothing 


Base-point 
positions 
in  x 


0.4 
.6 


1.0 


1.2 


1.4 
1.6 
1.8 


Scale 

variate   

z  FBI 

2.06  0.0048 

2.08  .0062 

2.29  .0149 

2.25  .0143 

2.32  .0151 

2.50  .0117 

2.50  .0117 

2.57  .0092 

2.83  .0010 

2.87  .0005 

3.00  .0000 

2.12 

2.13 

2.30 

2.30 

2.43 

2.49 

2.58 

2.75 

2.87 

2.03 
2.22 
2.30 
2.33 
2.43 
2.55 
2.57 
2.78 
2.80 
2.93 


From  table  A-3. 


Values  of  coefficients  1 
for  base  points 


YB2 


YB3 


YBA 


YB5 


-0.0373 
-.0483 
-.1144 
-.1090 
-.1163 
-.0977 
-.0977 
-.0813 
-.0202 
-.0139 
.0000 
.0089 
.0095 
.0150 
.0150 
.0138 
.0121 
.0088 
.0028 
.0005 


KB6 


0.9949 

0  049Q 

-U.UU04 

0.0001 

.9907 

.0588 

AA7C 

.0001 

.8566 

.0040 

.8948 

233? 

.0028 

.8245 

.3245 

.0050 

.5859 

.5859 

HQ77 

.0117 

.5859 

.5859 

0Q77 

.0117 

.4826 

.6858 

1 1  m 

.0138 

.1421 

.9536 

-  .yjouL, 

.0115 

.1026 

.9739 

-.0726 

.0095 

.0000 

1.0000 

.0000 

.0000 

-.0681 

.9780 

■  Woo 

-.U1Z5 

-.0726 

.9739 

1  f)9£! 

(lion 

-.1152 

.8462 

.6v  1 0 

-.U4  to 

-.1152 

.8462 

9Q7^ 

-.04  (O 

-.1101 

.6858 

-.0997 

.6005 

^71  Q 

-.0955 

-.0787 

.4678 

.oyyo 

-.1114 

-.0360 

.2332 

-.1090 

-.0139 

.1026 

Q7"3Q 

■  V  IOV 

-.0726 

.0025 

QQDQ 

.yyao 

.0207 

.0135 

-.1028 

.9197 

.1972 

.0150 

-.1152 

.8462 

.2975 

.0151 

-.1166 

.8132 

.3383 

.0138 

-.1101 

.6858 

.4826 

.0099 

-.0862 

.5122 

.6579 

.0092 

-.0813 

.4826 

.6858 

.0020 

-.0296 

.1972 

.9197 

.0016 

-.0256 

.1744 

.9344 

.0001 

-.0064 

.0508 

.9930 

x2 


Xa 


*5 


*6 


 Conventional  regression  symbols 


NOTES  FOR  PROGRAMS  B-l  AND  B-2 


YB1 


0.0004 
.0005 
.0043 
.0043 
.0092 
.0114 
.0140 
.0143 
.0095 

-.0026 
-.0296 
-.0478 
-.0554 
-.0813 
-.1071 
-.1101 
-.1028 
-.0976 
-.0429 


YB8  Y_ 

1.50 
1.69 
2.02 
2.20 
2.50 
2.50 
2.70 
2.89 
3.25 
3.58 
3.62 
3.66 
3.85 
3.55 
4.04 
3.88 

  3.60 

3.98 

  3.70 

  3.87 

0.0000  3.58 

.0020  3.69 

.0043  3.29 

.0054  3.56 

•0092  3.27 

.0133  3.10 

•0138  3.23 

.0135  2.74 

•  0128  2.98 

.0055  2.55 


X7 


x8 


X(I) 


SSX(I) 


Input  conventional  re- 
gression data  from 

^  tables  B-2  and  B-3. 

Sums  of  squares  and 
products. 


AA(I) 

MINV(AA,N,D,L,M) 


Matrix  to  be  inverted. 
System   library  routine 

for  matrix  inversion. 
Regression  coefficients' 
or  base  points. 
INTFOR(V,¥,YA,MS)  See  program  A-l 
INTSIX(V,Y,YP,MS)    See  program  A-3.' 


B(I] 
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PROGRAM  B-l.  —  Four-point  smoothing 


DIMENSION    X(30,10)  ,SSX(10,10)  , L  ( 10 ) , M < 1 0 )  , B ( 1 0) , Y I  NT ( 6 0 )  , Y AI NT (6 0 ) 
1, AA(6,6) 
COMMON/XX/  TAB(IGO) 

READ (5, 10 00) (<X(I,J) ,J=1,7) ,1=1,30) 
1000   FORMAT  ( 7F5 

DO  1002  I =1,10 

DO  1002  J=l,10 
1002  SSX (I, J)=0 .0 

DO   1003  1=1,7 

DO  1CQ4  J  =  I,7 

DO   1005  K=l,30 
100  5  SSX (I , J )  =  SSX ( I , J) +X ( K, I) *X ( K, J) 
1004  CONTINUE 
luO  3  CONTINUE 

DO  1006  1=1,7 

DO  1006  J=I,7 

1006  SSX ( J , I )  =  SSX ( I , J) 

WRITE (6,1010) ((SSX (I, J) ,1=1,7) ,J=1,7) 

1010  F ORM AT ( *  *,(7F12.8)) 
DO  1020  1=1,6 

DO  1020  J=l,6 
1020  A  A (I, J) -SSX (If  J) 

CALL  MINV ( AA,6,D,L  ,M) 

WRITE (5, 10 11)  ((AA(I,J)  ,1  =  1,6)  , J  =  l,6) 

1011  FORMAT (*  *,(6F13.8)) 
DO  1007  1=1,6 

B (I )  =  0  •  0 
DO  1003  J=l,6 
1008  B (I)=B(I) +AA (If  J) *SSX  (J, 7) 

1007  CONTINUE 

WRITE (6,10 11)  (B( I)  ,1  =  1,6) 
TAB(l) =0, k 
TAB (2)  =  0.  2 
DO  1012  1=3,8 

1012  TAB(I)=3( 1-2) 

V  =  0.  6 
MS=-1 

DO  1013  1=1,31 

CALL  INTFOR(Vf Yf YA,MS) 

Y INT(I) =Y 

YAINT(I)=YA 

\/  =  V  +  0.  T2 

1013  CONTINUE 

Y  AINT ( 1)  =  0  .  0 

WRITE (6,1014) (I ,YINT (I) ,1=1,31) 

1014  FORMAT (*   INTERPOLATED  FUNCTION*/ (8 ( 15, F10 . 6) )  ) 
WRITE(6,1015)  (I,YAINT(I)  ,1  =  1,31) 

1015  FORMAT ( *   INTERPOLATED  INTEGRAL*/ (8 (15, FIG . 6) ) ) 
STOP 

END 
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PROGRAM  B-2.  — Six-point  smoothing 


DIMENSION    X(30,1G) ,SSX(iO,10) , L(10) ,M(iO) ,B(1Q) , Y INT (60) , YPINT (60 ) 

_i,AA(8,S)   

COMMON /XX/  TAB (10  0 ) 

 READ  (5 » IQ  0  0)  ( (X  (I ,  J)  ,  J=l.,9)  ,1  =  1,30  )  

100  0   FORMAT  (9F5.it) 

00  1002  1  =  1,10  _____ 

00   100  2  J  =  1,10  ~    ~  -— —  — 

1Q  0  2  S  S  X  ( I ,  J )  =  0  .  0   

""DO   10  0  3  1  =1  ,9 

_00   100^  J=I,9   

DO   1005  K=l ,30 
10  05  SSX(I,J)=SSX(I,J) +X<K,I)*X<K, J) 
100  ^  CONTINUE  . 
100  3  CONTINUE  

DO   1006  1=1,9 

DO   1006  J=I,9 

1006  SSX  (J, I )  =  SSX (I, J) 

WRITE (6, 10 10) ((SSX(I,J),I=1,9) f J=lf9__ 
101 0   FORMAT (*  »,(9F12.8>) 

DO   1020  1=1,8   

DO   1023  J=l,8 

1020   AA ( I , J) =SSX (I, J)   

CALL   MINV ( AA,8,D,L,M) 

WRITE (6 ,1 011) ( ( AA ( I , J) ,1=1 , 8) ,J=i,8) 
lOlT  FORMAT (*  *,(8F13.8)) 

DO      10)7  1=1,8 

D  ( I )  =  0  •  0 

DO   1003  J=l,8 
1  0  0  8  9  (I)=B(I)  +AA(I,J)*SSX(J,9) 

1007  CONTINUE 

WRITE (6 ,1011) (B(I),I=1,8) 

TAB  (1)  =  0. k 

TAB  (2)  =  0 . 2 

DO  1012  1=3,10 

1012  tap (i) =3( i-2) 

V=C  .8 
MS=-1 

DO   1013  1=1,31 

CALL   INTSIX(V,Y,YP ,MS) 

Y INT (I ) =Y 

YPINT(I)=YP 

V=V  +  C . 02 

1013  CONTINUE 

WRITE i6,101h)  (I, Y I  NT  (I) ,1  =  1,31) 
101^  FORMAT  ( *    INTERPOLATED   FUNCTION*/ ( 8 ( 15 , F 10 . 6 ) ) ) 

WRITE (6, 10 15) (I, YPINT (I) ,1=1,31) 
1015  cORMAT  (*      INTERPOLATED  DER I VA T I VE* / ( 8 ( I  5 »  Fi 0  . 6) ) ) 

STOP 

END 


*  1976  GP0  1750-S-671 -583/49 
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